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ABSTRACT 

We report the discovery of a massive {Mp = 9.04 ± 0.50 Mjup) planet transiting the bright {V = 8.7) 
F8 star HD 147506, with an orbital period of 5.63341±0.00013 days and an eccentricity of e = 0.520± 
0.010. From the transit light curve we determine that the radius of the planet is Rp = 0.982±q ]'o5 i?jup- 
HD 147506b (also coined HAT-P-2b) has a mass about 9 times the average mass of previously-known 
transiting exoplanets, and a density of pp = 11.9 gcm""^, greater than that of rocky planets like 
the Earth. Its mass and radius are marginally consistent with theories of structure of massive giant 
planets composed of pure H and He, and may require a large (> lOOMgg) core to account for. The high 
eccentricity causes a 9-fold variation of insolation of the planet between peri- and apastron. Using 
follow-up photometry, we find that the center of transit is Tmid = 2,454,212.8559 ± 0.0007 (HJD), and 
the transit duration is 0.177 ± 0.002 d. 

Subject headings: stars: individual: HD 147506 - planetary systems: individual: HD 147506b, HAT- 
P-2b 



1. INTRODUCTION 

To date 18 extrasolar planets have been found which 
transit their parent stars and thus yield values for their 
mass and radius^'^. Masses range from 0.3 Mj to about 
1.9 Mj, and radii from 0.7i?j to about 1.4i?j. The ma- 
jority fit approximately what one expects from theory 
for irradiated gas giant planets (e.g. Fortney, Marley, 
& Barnes 2006, and references therein), although there 
are exceptions: HD 149026b has a small radius for its 
mass (Sato et al. 2005), implying that it has a large 
heavy core (~70Af®; Laughlin et al. 2005), and sev- 
eral (HD 209458b, HAT-P-lb, WASP-1) have unexpect- 
edly large radii for their masses, perhaps suggesting some 
presently unknown source of extra internal heating (Guil- 
lot & Showman 2002; Bodenheimer, Laughlin, & Lin 
2003). The longest period and lowest density transit- 
ing exoplanet (TEP) detected so far is HAT-P-lb with 
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P = 4.46 d (Bakos et al. 2007). AU TEPs have orbits 
consistent with circular Keplerian motion. 

From existing radial velocity (RV) data, it might be 
expected that there are some close-in (semi-major axis 
< 0.07 AU, or P < 10 d) giant planets with masses con- 
siderably larger than any of the 18 transiting planets 
now known. A well-known example, considering only ob- 
jects below the Deuterium burning threshold (~ 13 Mj, 
e.g. Burrows et al. 1997), is r Boo b, which was de- 
tected from RV variations, and has a minimum mass of 
Mp sin 1=3.9 Mj and orbits only 0.046 AU from its star 
(Butler et al. 1997). Another example is HIP14810 b 
(Wright et al. 2007) with similar mass, and orbital pe- 
riod of 6.7 d and semi-major axis of 0.069 AU. At this 
orbital distance the a priori probability of such a planet 
transiting its star is about 10%. Thus "super-massive" 
planets should sometimes be found transiting their par- 
ent stars. We report here the detection of the first such 
TEP, and our determination of its mass and radius. This 
is also the longest period TEP, and the first one to ex- 
hibit highly eccentric orbit. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Detection of the transit in the HATNet data 

HD 147506 is an F8 star with visual magnitude 8.7 
and Hipparcos parallax 7.39 ± 0.88 mas (Ferryman et al. 
1997). It was initially identified as a transit candidate in 
our internally labeled field G193 in the data obtained by 
HATNet'si"* (Bakos et al. 2002, 2004) HAT-6 telescope 
at the Fred Lawrence Whipple Observatory (FLWO) of 
the Smithsonian Astrophysical Observatory (SAO). The 
detection of a '--^5mmag transit with a 5.63 d period in 
the light curve consisting of ^7000 data-points (with a 
5.5 min cadence) was marginal. Fortunately the star was 
in the overlapping corner with another field (G192) that 
has been jointly observed by HATNet's HAT-9 telescope 
at the Submillimeter Array (SMA) site atop Mauna Kea, 

http://www.hatnet.hu 
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Fig. 1. — The upper panel shows the unbinned HATNet and WHAT joint hght curve with 26400 data-points, phased with the P = 5.63341 d 
period. The 5 mmag deep transit is detected with a signal-to-noise of 26. The middle panel shows the same HATNet and WHAT data 
zooming in on the transit and binned with a </> = 0.0005 bin-size. The lower panel displays the Sloan z-band photometry taken with the 
FLWO 1.2 m telescope. Over-plotted is our best fit obtained with the Mandcl & Agol (2002) formalism. 



Hawaii, and by the Wise HAT telescope (WHAT, Wise- 
Observatory, Israel; Shporer et al. 2006), for an extended 
period that yielded ^6700 and ^3900 additional data- 
points, respectively. The transit was independently de- 
tected and confirmed with these data-sets. By chance 
the candidate is in yet another joint field (G149) of HAT- 
Net (HAT-7 at FLWO) and WHAT, contributing ~6200 
and ~2200 additional data points, respectively. Alto- 
gether this resulted in a light curve with exceptional 
time-coverage (570 days), an unprecedented number of 
data-points (26461 measurements at 5.5 min cadence), 
and an rms of 5 mmag. It is noteworthy that the network 
coverage by WHAT (longitude 35 E), HATNet at FLWO 
(111 W) and HATNet at Hawaii (155 W) played an im- 
portant role in detecting such a long period and shallow 
transit. Data were reduced using astrometry from Pal 
& Bakos (2006), and with a highly fine-tuned aperture 
photometry. We applied our external parameter decor- 
relation (EPD) technique on the light curves, whereby 
deviations from the median were cross-correlated with a 
number of "external parameters" , such as the X and Y 
sub-pixel position, hour-angle, and zenith distance. We 
have also applied the Trend Filtering Algorithm (TFA; 
Kovacs, Bakos, & Noyes 2005, hereafter KBN05) along 
with the Box Least Squares (BLS: Kovacs, Zucker, & 
Mazeh 2002) transit-search algorithm in our analysis. 
TFA and BLS were combined in signal-reconstruction 
mode, assuming general signal shape, as described in 
KBN05. The detection of this relatively shallow transit 
is a good demonstration of the strengths of TFA. The 
upper panel of Fig. 1 shows the unbinned light curve with 



all 26400 data points, whereas the middle panel displays 
the transit binned to 1/2000 of the period (4 minutes). 

We note that due to the large amount of data, the binned 
light curve is of similar precision as a single-transit ob- 
servation by a Im-class telescope. 

After several failed attempts (due to bad weather and 
instrumental failure) to carry out high-precision photo- 
metric follow-up observations from FLWO, Wise Ob- 
servatory, Konkoly Observatory, and the Clay Center 
(Boston), we finally succeeded in observing a full tran- 
sit using the Kepler Cam detector on the FLWO 1.2 m 
telescope (see Holman et al. 2007) on UT 2007 April 
22. The Sloan z-band light curve is shown in the lower 
panel of Fig. 1. From the combined HATNet and Ke- 
plerCam photometry, spanning a baseline of 839 days, 
we derive a period of 5.63341 ± 0.00013 d and an epoch 
of mid-transit of T^d = 2,454,212.8559±0.0007d (HJD). 
From the FLWO 1.2 m data alone (and the analytic light 
curve fit as described later), the length of transit is 
0.177±0.002d (4 hours, 15 minutes), the length of ingress 
is 0.012 ± 0.002 d (17.5 minutes), and the depth (at the 
middle of the transit) is 0.0052 mag. 

2.2. Early spectroscopy follow-up 

Initial follow-up observations were made with the CfA 
Digital Speedometer (DS; Latham 1992) in order to char- 
acterize the host star and to reject obvious astrophysi- 
cal false-positive scenarios that mimic planetary tran- 
sits. These observations yielded values of T^s = 6250 K, 
log 5 = 4.0 and osini = 22kms~^, corresponding to a 
moderately-rotating main sequence F star. The radial 
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TABLE 1 
Radial Velocities for HD 147506. 



TABLE 2 

Summary of stellar parameters for HD 147506. 
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Parameter 



Value 



Source 



^The RVs include the barycentric correction. 

^Only the Keck and Lick data-points are shown here. Consult the 
electronic edition for a full data-set that includes the CfA DS measure- 
ments. 

velocity (RV) measurements showed an rms residual of 
~0.82kms~^, slightly larger than the nominal DS pre- 
cision for a star with this rotation, and suggested that 
they may be variable. With a few dozen additional DS 
observations, it was found that the RV appeared peri- 
odic with P » 5.63d, semi-amplitude ~lkms^^, and 
phasing in agreement with predictions from the HAT- 
Net+WHAT light curve. This gave strong evidence that 
there really was an RV signal resulting from Keplerian 
motion, although the precision was insufficient to estab- 
lish the orbit with confidence. Altogether we collected 
53 individual spectra spanning a time-base of more than 
a year (Table 1). 

2.3. High-precision spectroscopy follow-up 

In order to confirm or refute the planetary nature of 
the transiting object, we pursued follow-up observations 
with the HIRES instrument (Vogt et al. 1994) on the 
W. M. Keck telescope and with the Hamilton Echelle 
spectrograph at the Lick Observatory (Vogt 1987). The 
spectrometer slit used at Keck is (y.'86, yielding a re- 
solving power of about 55,000 with a spectral coverage 
between about 3200 and 8800 A. The Hamihon Echelle 
spectrograph at Lick has a similar resolution of about 
50,000. These spectra wcirc; used to i) more fully charac- 
terize the stellar properties of the system, ii) to obtain a 
radial velocity orbit, and to iii) check for spectral line bi- 
sector variations that may be indicative of a blend. We 
gathered 13 spectra at Keck (plus an iodine- free tem- 
plate) spanning 207 days, and 10 spectra at Lick (plus 
template) spanning 50 days. The radial velocities mea- 
sured from these spectra are shown in Table 1, along with 
those from the CfA DS. 



logg 

i;sini(kms~'^) 
[Fe/H](dex) 
Distance (pc) 
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Mass (Mq) 
Radius (Rq) 
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A spectral synthesis modeling of the iodine-free Keck 
template spectrum was carried out using the SME soft- 
ware (Valenti & Piskunov 1996), with the wavelength 
ranges and atomic line data described by Valenti & Fis- 
cher (2005). Results are shown in Table 2. The val- 
ues obtained for the effective temperature (Tch), surface 
gravity (logf/), and projected rotational velocity (wsini) 
are consistent with those found from the CfA DS spectra. 
As a check on Teg, we collected all available photometry 
for HD 147506 in the Johnson, Cousins, 2MASS, and Ty- 
cho systems, and applied a number of color- temperature 
calibrations (Ramirez & Melendez 2005; Masana et al. 
2006; Casagrande et al. 2006) using 7 different color 
indices. These resulted in an average temperature of 
'~6400±100K, somewhat higher than the spectroscopic 
value but consistent within the errors. 

Based on the Hippareos parallax (tt = 7.39±0.88mas), 
the apparent magnitude V = 8.71 ± 0.01 (Droege et al. 
2006), the SME temperature, and a bolometrie correc- 
tion of BCv = -0.011 ± 0.011 mag (Flower 1996), appli- 
cation of the Stefan-Boltzmann law yields a stellar radius 
of i?* = 1.84 ± 0.24 i?Q. 

A more sophisticated approach to determine the stellar 
parameters uses stellar evolution models along with the 
observational constraints from spectroscopy. For this we 
used the models by Yi et al. (2001) and Demarque et 
al. (2004), and explored a wide range of ages to find all 
models consistent with Tcs, My, and [Fe/H] within the 
observational errors. Here My — 3.05 ± 0.26 is the abso- 
lute visual magnitude, as calculated from V and the Hip- 
pareos parallax. This resulted in a mass and radius for 

f0 and i?„ = 1.85 ±[^:iJ i?0, 



the star of = 1.42±;]:}§ Mq and i?„ = 1.85±[^:^J Rr, 
and a best-fit age of 2.7lQ;g Gyr. Other methods that rely 
on the Hippareos parallax, such as the Padova^^ stellar 
model grids (Girardi et al. 2002), consistently yielded a 
stellar mass of ~ 1.4 and stellar radius ^ 1.8 i?0. 

If we do not rely on the Hippareos parallax, and use 
log 5 as a proxy for luminosity (instead of My), then the 
Y^ stellar evolution models yield a smaller stellar mass of 
= 1.29±E;;1^ Mq, radius of i?„ = 1.46 ±0:^6 R^^ and 
best fit age of 2.6l2'5 ^Y^- The surface gravity is a sen- 
sitive measure of the degree of evolution of the star, as is 
luminosity, and therefore has a very strong influence on 
the radius. However, logg is a notoriously diflacult quan- 
tity to measure spectroscopically and is often strongly 
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correlated with other spectroscopic parameters. 

It has been pointed out by Sozzctti et aL (2007) that 
the normaUzed separation a/i?* can provide a much bet- 
ter constraint for stellar parameter determination than 
log 5. The a/Ri, quantity can be determined directly 
from the photometric observations, without additional 
assumptions, and it is related to the density of the cen- 
tral star. As discussed later in § 6, an analytic fit to 
the FLWO 1.2 m light curve, taking into account an ec- 
centric orbit, yielded a/i?,* = S.TT^qqj. Using this as 
a constraint, along with Toff and [Fe/H], we obtained 

0.06 A/f.. r? - 1 /17 -1-0.04 



= 1.30 ±l-^il Mq 



= 1.47 ±[j:^^ Rq and age of 



2.6;?:^ Gyr. The logg = 4.214±g:gf| derived this way is 
consistent with former value from SME. 

As seen from the above discussion, there is an inconsis- 
tency between stellar parameters depending on whether 
the Hipparcos parallax is employed or not. Methods rely- 
ing on the parallax (Stefan-Boltzmann law, stellar evolu- 
tion models with My constraint, etc) prefer a larger mass 
and radius 1.4^^0, ^ 1.8 i?0, respectively), whereas 
methods that do not rely on the parallax (stellar evo- 
lution models with logg or a/i?* constraint) point to 
smaller mass and radius (~ 1.3 Mq, ~ 1.46 i?0, respec- 
tively). We have chosen to rely on the a/R^ method, 
which yields considerably smaller uncertainties and a cal- 
culated transit duration that matches the observations. 
Additionally, it implies an angular diameter for the star 
(0 = 0.127^0.014 Bias) that is in agreement with the 
more direct estimate of ^ = 0.117 ± 0.001 mas from 
the near-infrared surface-brightness relation by Kervella 
et al. (2004). The later estimate depends only on the 
measured V — Ks color and apparent Kg magnitude (ig- 
noring extinction) from 2MASS (Skrutskie et al. 2006), 
properly converted to the homogenized Bessell & Brett 
system for this application (following Carpenter 2001). 
We note that our results from the a/i?^ method imply a 
somewhat smaller distance to HD 147506 than the one 
based on the Hipparcos parallax. The final adopted stel- 
lar parameters are listed in Table 2. 

3.1. Stellar jitter 

Stars with significant rotation are known to exhibit 
excess scatter ("jitter") in their radial velocities (e.g., 
Wright 2005, and references therein), due to enhanced 
chromospheric activity and the associated surface inho- 
mogcneitics (spottedncss). This jitter is in addition to 
the internal errors in the measured velocities, and could 
potentially be significant in our case. We note that af- 
ter pre-whitening the light curve with the transit com- 
ponent, we found no significant sinusoidal signal above 
0.3mmag amplitude. From this we conclude that there 
is no very significant spot activity on the star (in the 
observed 500 day window). In order to estimate the 
level of chromospheric activity in the star, we have de- 
rived an activity index from the Ca II H and K lines 
in our Keck spectra of logi?^^^ = —4.72 ± 0.05. For 
this value the calibration by Wright (2005) predicts ve- 
locity jitter ranging from 8 to 16 ms^^. An earlier cali- 
bration by Saar, Butler, & Marcy (1998), parametrized 
in terms of the projected rotational velocity, predicts a 
jitter level of up to 50ms~"'^ for our measured vsini of 
20kms~^. A different calibration by the same authors in 
terms of R'^k gives 20ms~^. An additional way to esti- 
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Fk;. 2. — The upper panel shows the RV measurements phased 
with the period of P = 5.63341 d. The zero-point in phase cor- 
responds to the epoch of mid- transit. Large filled circles indicate 
Keck and Lick points. Small open circles denote CfA DS data (not 
used for the fit). Overlaid is the fit that was based only on the 
Keck and Lick data assuming 60ms~^ stellar jitter. The middle 
panel shows the residuals from the fit. The lower panel exhibits 
the line bisector spans on the same scale as the upper panel. No 
variation in the line bisectors is seen concomitant with that in the 
RVs, essentially confirming the planetary nature of the transiting 
object. 

mate the jitter is to compare HD 147506 to stars of the 
Lick Planet Search program (Gumming, Marcy, & Butler 

1999) that have similar properties (0.4 < B — V < 0.5, 
vsini > 15kms^^). There are four such stars (J. John- 
son, private communication), and their average jitter is 
45ms~^. A more direct measure for the particular case 
of HD 147506 may be obtained from the multiple ex- 
posures we collected during a 3-night Keck run in 2007 
March. Ignoring the small velocity variations due to or- 
bital motion during any given night, the overall scat- 
ter of these 8 exposures relative to the nightly means is 
~20ms~^. This may be taken as an estimate of the jit- 
ter on short timescales, although it could be somewhat 
larger over the entire span of our observations. Alto- 
gether, it is reasonable to expect the jitter to be at least 
10ms~^, and possibly around 30 50 ms^^ for this star. 

4. SPECTROSCOPIC ORBITAL SOLUTION 

We have three velocity data sets available for analysis: 
13 relative radial velocity measurements from Keck, 10 
from Lick, and 53 measurements from the CfA DS, which 
are nominally on an absolute scale (Table 1). Given the 
potential effect of stellar jitter, we performed weighted 
Keplerian orbital solutions for a range of jitter values 
from 10 to 80ms~^ with 10ms~^ steps. These jitter 
values were added in quadrature to all individual inter- 
nal errors. We performed separate fits for the star or- 
bited by a single planet, both with and without the CfA 
DS measurements, since these have errors (^600 ms^^) 
significantly larger than Keck (5 9nis^^) or Lick (40 
90ms~^). In all of these solutions we held the period 
and transit epoch fixed at the photometric values given 
earlier. The parameters adjusted are the velocity semi- 
amplitude K, the eccentricity e, the longitude of perias- 
tron u), the center-of-mass velocity for the Keck relative 
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velocities 7, and offsets IS.vkl between Keck and Lick 
and /\vKC between Keck and CfA DS. The fitted pa- 
rameters were found to be fairly insensitive to the level 
of jitter assumed. However, only for a jitter of ^GOms"^ 
(or ~70ms~^ when the CfA DS data are included) did 
the approach values expected from the number of de- 
grees of freedom. There arc thus two possible conclu- 
sions: if we accept that HD 147506 has stellar jitter at 
the 60 ms^^ level, then a single-planet solution such as 
ours adequately describes our observations. If, on the 
other hand, the true jitter is much smaller (< 20ms~^), 
then the extra scatter requires further explanation (see 
below) . Our adopted orbital parameters for the simplest 
single-planet Keplerian solution are based only on the 
more precise Keck and Lick data, and assume the jitter 
is 60ms~^ (Table 3). The orbital fit is shown graphically 
in the upper panel of Fig. 2. In this figure, the zero-point 
of phase is chosen to occur at the epoch of mid-transit, 
Trnid = 2,454,212.8559 (HJD). The most significant re- 
sults are the large eccentricity (e = 0.520 ± 0.010), and 
the large velocity semi-amplitude {K = 1011 ± 38 ms~^, 
indicating a very massive companion). As we show 
in the next section (§ 5), the companion is a planet, 
i.e. HD 147506b, which we hereafter refer to as HAT- 
P-2b. 

As a consistency check we also fitted the orbits by fix- 
ing only the period, and leaving the transit epoch as a 
free parameter. We found that for all values of the stellar 
jitter the predicted time of transit as derived from the RV 
fit was consistent with the photometric ephemeris within 
the uncertainties. We also found that in these fits the 
orbital parameters were insensitive to the level of jitter 
and to whether or not the CfA DS data were included. 
The eccentricity values ranged from 0.51 to 0.53. 

4.1. Solutions involving two planets 

If we assimie that the true stellar jitter is small, then 
the excess scatter in the RV fit could be explained by a 
third body in the system, i.e., a hypothetical HAT-P-2c. 
In addition, such a body could provide a natural dynam- 
ical explanation for the large eccentricity of HAT-P-2b at 
this relatively short period orbit. Preliminary two- planet 
orbital fits using all the data yielded solutions only sig- 
nificant at the 2-sigma level, not compelling enough to 
consider as evidence for such a configuration. Additional 
RV measurements are needed to firmly establish or refute 
the existence of HAT-P-2c. 

We also exploited the fact that the HATNet light curve 
has a unique time coverage and precision, and searched 
for signs of a second transit that might be due to an- 
other orbiting body around the host star. Successive 
box-prewhitening based on the BLS spectrum and as- 
suming trapezoidal-shape transits revealed no secondary 
transit deeper than the 0.1% level and period < 10 days. 

5. EXCLUDING BLEND SCENARIOS 

As an initial test to explore the possibility that the 
photometric signal we detect is a false positive (blend) 
due to c;ontamination from an unresolved eclipsing bi- 
nary, we modeled the light curve assuming there are three 
coeval stars in the system, as described by Torres et al. 
(2004). We were indeed able to reproduce the observed 
light curve with a configuration in which the brighter ob- 
ject is accompanied by a slightly smaller F star which is 



in turn being eclipsed by a late-type M dwarf. However, 
the predicted relative brightness of the two brighter ob- 
jects at optical wavelengths would be ~ 0.58, and this 
would have been easily detected in our spectra. This 
configuration can thus be ruled out. 

The reality of the velocity variations was tested by 
carefully examining the spectral line bisectors of the 
star in our more numerous Keck spectra. If the veloc- 
ity changes measured are due only to distortions in the 
line profiles arising from contamination of the spectrum 
by the presence of a binary with a period of 5.63 days, 
we would expect the bisector spans (which measure line 
asymmetry) to vary with this period and with an am- 
plitude similar to the velocities (see, e.g., Queloz et al. 
2001; Torres et al. 2005). The bisector spans were com- 
puted from the cross-correlation function averaged over 
15 spectral orders blucward of 5000 A and miaffected by 
the iodine lines, which is representative of the average 
spectral line profile of the star. The cross correlations 
were performed against a synthetic spectrum matching 
the effective temperature, surface gravity, and rotational 
broadening of the star as determined from the SME anal- 
ysis. As shown in Fig. 2, while the measured velocities 
exhibit significant variation as a function of phase (up- 
per panel), the bisector spans are essentially constant 
within the errors (lower panel). Therefore, this analysis 
rules out a blend scenario, and confirms that the orbiting 
body is indeed a planet. 

6. PLANETARY PARAMETERS 

For a precise determination of the physical properties 
of HAT-P-2b we have modeled the FLWO 1.2 m Sloan 
2r-band photometric data shown in Fig. 1. The model is 
an eccentric Keplerian orbit of a star and planet, thus 
accounting for the nonuniform speed of the planet and 
the reflex motion of the star. Outside of transits, the 
model flux is unity. During transits, the model flux is 
computed using the formalism of Mandel & Agol (2002), 
which provides an analytic approximation of the flux of 
a limb-darkened star that is being eclipsed. The free pa- 
rameters were the mid-transit time Tmid, the radius ratio 
Rp/ Ri,, the orbital inclination i, and the scale parameter 
a/Ri,, where a is the semimajor axis of the relative orbit. 
The latter parameter is determined by the time scales of 
the transit (the total duration and the partial-transit du- 
ration), and is related to the mean density of the star (see 
§ 3). The orbital period, eccentricity, and argument of 
pericenter were flxed at the values determined previously 
by fitting the radial velocity data. The limb darkening 
law was assumed to be quadratic, with coefficients taken 
from Claret (2004). 

To solve for the parameters and their imccrtainties, we 
used a Markov Chain Monte Carlo algorithm that has 
been used extensively for modeling other transits (see, 
e.g. Winn, Holman, & Fuentes 2007; Holman et al. 2007). 
This algorithm determines the a posteriori probability 
distribution for each parameter, assuming independent 
( "white" ) Gaussian noise in the photometric data. How- 
ever, we found that there are indeed correlated errors. 
Following Gillon et al. (2006), we estimated the red noise 
(Jr via the equation 

= l-l/AT ' 
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TABLE 3 



Orbital fit and planetary parameters for the HAT-P-2 

SYSTEM. 


Parameter 


Value 


Period (d)^ 
T„M (HJD)=» 
Transit duration (day) 
Ingress duration (day) 


5.63341 ± 0.00013 
2,454,212.8559 ± 0.0007 
0.177 ±0.002 
0.012 ± 0.002 


Stellar jitter (ms^^)" 
7 (ms"-'-)'^ 
K (ms-i) 
(deg) 

e 

Tperi (HJD) 

Avkl (ms-i) 
AvKC (kms-i)d 
/(M) (Mq) 
Mpsini (Mj) 
a* sini (km) 


60 

-278 ± 20 
1011 ±38 
179.3 ±3.6 
0.520 ±0.010 
2,454,213.369 ±0.041 
—380 ± 35 
19.827 ±0.087 
(376 ±42) X 10-9 
7.56 ± 0.28([M* + Mp]/M0])2/3 
(0.0669 ± 0.0025) x 10*^ 


a,ei (AU) 

ip(deg) 

Mp{Mj) 

ppigcm ^) 


0.0677 ±0.0014 
> 84.6° (95% confidence) 
9.04 ±0.50 
0.982±!]0g8 
11.9±tj 
227±f| 


'^Fixed in the orbital fit. 
^Adopted (sec text). 

'^The 7 velocity is not in an absolute reference frame. 
"^The offset between Keck and CfA DS is given for reference from a 
fit that includes all data sets, but does not affect our solution. 



where txi is the standard deviation of the out-of-transit 
flux of the original (unbinned) hght curve, cttv is the stan- 
dard deviation of the hght curve after binning into groups 
of N data points, and iV = 40 corresponds to a bin- 
ning duration of 20 minutes, which is the ingress/egress 
time scale that is critical for parameter estimation. With 
white noise only, = axj^pN and Ur = 0. We added Or 
in quadrature to the error bar of each point, effectively 
inflating the error bars by a factor of 1.25. 

The result for the radius ratio is Rp/Ri, = 0.0684 ± 
0.0009, and for the scale parameter a/Ri, — 9.77tQQ2- 
The a posteriori distribution for a/ R^, is very asymmet- 
ric because the transit is consistent with being equatorial: 
i > 84?6 with 95% confldence. We conflrmed that these 
uncertainties are dominated by the photometric errors, 
rather than by the covariances with the orbital param- 
eters e, iv, and P, and hence we were justifled in flx- 
ing those orbital parameters at constant values. Based 
on the inclination, the mass of the star (Table 2) and 
the orbital parameters (Table 3), the planet mass is then 
9.04 ± 0.50 Mj. Based on the radius of the star (Table 2) 
and the above Rp/R^, determination, the radius of the 
planet is Rp = 0.982 ip;?!]! These properties are 

summarized in Table 3. 



7. DISCUSSION 

In comparison with the other 18 previously known 
transiting exoplanets, HAT-P-2b is quite remarkable 
(Fig. 3, Fig. 4). Its mass of 9.04 ± 0.50 Mj is ~5 times 
greater than any of these 18 other exoplanets. Its mean 
density p = 11.9±f g gcm"'^ is ~9 times that of the dens- 
est known exoplanet (0GLE-TR-113b, p = 1.35 gcm"^) 
and indeed greater than that of the rocky planets of the 
Solar System {p — 5.5 gcm"'^). Its surface gravity of 



0.1 1 10 100 

Mp [Mjyp] 

Fig. 3. — The mass-radius diagram of known TEPs (from 
www.exoplanet.eu and references therein), Jupiter and Saturn 
(large filled circles), and low-mass stars from Beatty et al. (2007). 
HAT-P-2b is an intermediate mass object that is still in the plan- 
etary regime (well below 13Afj). Overlaid are equidensity lines 
(labeled), Baraffe et al. (1998) (stellar) and Barafll'e et al. (2003) 
(zero insolation planetary) isochrones for ages of 0.5 Gyr (upper, 
dotted line) and 5 Gyr (lower dashed-dotted line) , respectively. 

227 ±f| ms^^ is 7 times that of any of the previously 
known TEPs, and 30 times that of HAT-P-lb (Fig. 4). 

We may compare the mass and radius for HAT-P-2b 
with evolutionary models, including irradiation, as re- 
cently presented by Fortney, Marley, & Barnes (2006) 
(hereafter FMB06). Given the inferred stellar luminos- 
ity (Table 2), and the time-integral of the insolation over 
an entire period (taking into account the orbital parame- 
ters, notably e and Orei), the equivalent semi- major axis 
tteq for the same amount of irradiation if the central star 
were solar is 0.036 AU. At that separation, FMB06 find 
for a pure hydrogen/helium planet of mass 9Mj and 
age of 4.5 Gyr a planetary radius about 1.097 A 
100 core has a negligible effect on the radius (yield- 
ing 1.068 which is not surprising, since the mass of 
such a core is only a few percent of the total mass. For 
younger ages of 1 Gyr and 0.3 Gyr the radii are larger: 
1.159 i?j and 1.22 i?j for coreless models, respectively. 
Our observed radius of 0.982 i?j is smaller than any of 
the above values (4.5, 1, 0.3 Gyr, with or without 100 
core). Since the let positive error-bar on our radius deter- 
mination is 0.038 the inconsistency is only marginal. 
Nevertheless, the observed radius prefers either larger age 
or bigger core-size, or both. Given the age of the host 
star (2.6 Gyr, Table 2) the larger age is an unlikely expla- 
nation. The required core-size for this mass and radius 
according to FMB06 would be 300 , which amount of 
icy and rocky material may be hard to account for. 

Fig. 3 also shows a theoretical mass-radius relation for 
objects ranging from gas giant planets to stars (Baraffe 
et al. 1998, 2003). Note that HAT-P-2b falls on the re- 
lation connecting giant planets to brown dwarfs to stars. 
It thus appears to be intermediate in its properties be- 
tween Jupiter-like planets and more massive objects like 
brown dwarfs or even low mass stars. According to the- 
ories, stars with mass > 0.2 Mq have a core, where in- 
ternal pressure is dominated by classical gas (ions and 
electrons), and the R (x M radius-mass relation holds 
in hydrostatic equilibrium (for a review and details on 
the following relations see e.g. Chabrier et al. 2000). Be- 
low ~ 0.075 M0 (80 Mj) mass, however, the equation 
of state in the core becomes dominated by degenerate 
electron gas (R oc M^^^^ for full degeneracy), yield- 



A Super-Massive Transiting Planet 



7 



100 



E 







HAT-P-2b * 

- 


{ 


f * I 






f . ^ 

] 


f 



1 2 3 4 5 6 

P [days] 



Fig. 4. — The surface gravity of TEPs as a function of orbital 
period. Data taken from Southworth, Wheatlcy, & Sams (2007) 
with the exception of HAT-P-2b. This object is clearly not obeying 
the suspected correlation of gp and P for "Jupiter-mass" objects. 

ing an expected minimum in the mass-radius relation- 
ship (around 73 Mj). Below this mass, the partial de- 
generacy of the object and the classical [R oc M^/^) 
Coulomb pressure together yield an almost constant ra- 
dius {R cx M^i/*). HAT-P-2b is a demonstration of 
this well known phenomenon. (The approximate relation 
breaks below M ~ 4 Mj, where the degeneracy saturates, 
and a classical mass-radius behaviour is recovered). 

Compared to the other 18 known transiting planets, 
HAT-P-2b is also unique in having an orbit with remark- 
ably high eccentricity. The primary question is how such 
an eccentricity was created in the first place. One pos- 
sible explanation could be that the planet was scattered 
inward from a larger orbit, acquiring a high eccentricity 
in the process (Ford & Rasio 2007; Chatterjee, Ford, & 
Rasio 2007). If so, then the scattering event might have 
caused its new orbital plane to be inclined relative to the 
plane of the original disk, and hence out of the equatorial 
plane of the parent star (e.g. Fabrycky & Tremaine 2007). 
This angle between these two planes should be readily 
measurable from the Rossiter-McLauglilin effect (Winn 
et al. 2005). Indeed, the star HD 147506 is an ideal sub- 
ject for studying this effect, because its rapid rotation 
should lead to a relatively large Rossiter-McLaughlin sig- 
nal. 

There arc a number of other interesting issues related 
to the high eccentricity of HAT-P-2b. During its 5.63 day 
orbit, the insolation reaching the planet's surface varies 
by a factor of 9. Assuming an albedo of 0.1 (Rowe et al. 
2006) and complete redistribution of insolation energy 
over the surface of the planet, the equilibrium tempera- 
ture varies from about 2150 K at pcriastron to 1240 K at 
apastron. This would have a major influence on atmo- 
spheric dynamics and photochemistry. 

It is interesting to compare the properties of the HAT- 
P-2 system with the r Boo system, which - as already 
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noted - harbors a close-in planet with minimum mass 
Mp sin i = 3.9A'/j. Similarities of the two parent stars 
include the nearly identical masses, effective temper- 
ature, and the rapid rotation, although r Boo, with 
[Fe/H\ = -1-0.28, is somewhat more metal rich than 
HD 147506, with [Fe/H] = +0.12. A striking differ- 
ence is that, while the orbital eccentricity of HAT-P-2b 
is 0.5, the eccentricity of r Boo b is not measurably dif- 
ferent from zero. However, r Boo b's orbital period, 3.3 
days, is almost half that of HAT-P-2b. A large fraction 
of close-in planets with 5 < P < 10 days have signifi- 
cant eccentricities (0.1 < e < 0.3) although not as large 
as HAT-P-2b. For discussion on the eccentricity distri- 
bution see Juric & Tremaine (2007). As circularization 
timescales are thought to be very steep functions of the 
orbital semi-major axis (Terquem et al. 1998), one could 
then argue that HAT-P-2b's large value of e is due to 
either the fact that the planet's orbit is not yet circu- 
larized (while T Boo b's instead is), or to the presence 
of a second planet in the HAT-P-2 system, or to rather 
different formation/migration scenarios altogether. 

HD 147506, with visual magnitude 8.71, is the foiuth 
brightest among the known stars harboring transiting 
planets. Therefore it has special interest because of 
the possibilities for foUowup with large space or ground- 
based telescopes. 
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